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Abstract 

Loop quantum cosmology allows for arbitrary superpositions of the 
triad variable. We show here how these superpositions can become in- 
distinguishable from a classical mixture by the interaction with fermions. 



5-H . 

We calculate the reduced density matrix for a locally rotationally sym- 
metric Bianchi I model and show that the purity factor for the triads 
decreases by decoherence. In this way, the Universe assumes a definite 
orientation. 



^Address from October 1, 2012: Max Planck Institute for Gravitational Physics (Albert 
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1 Introduction 



Quantum theory seems to be a universal framework for all interactions. As 
such, it should be applicable to the Universe as a whole, leading to quantum 
cosmology. Since gravity is the dominant interaction on large scales, the 
formalism of quantum cosmology must be based on a theory of quantum 
gravity. 

A direct quantization of general relativity in the canonical formalism leads 
to quantum geometrodynamics and the Wheeler-DeWitt equation [1]. An 
alternative canonical approach is loop quantum gravity [H [2] , which leads to 
loop quantum cosmology when applied to cosmological models [3]. In loop 
quantum cosmology, the quantum state obeys a difference equation, and 
not a differential equation. The reason for this is the presence of discrete 
spectra for geometric operators (such as the area operator) in the full theory. 
This discreteness seems to facilitate the avoidance of the classical big-bang 
singularity. 

In loop quantum cosmology, one does not work with three-dimensional 
metrics, but with triads. This leads to an important difference to quantum 
geometrodynamics: triads can have two different orientations (left-handed 
and right-handed), which both lead to the same metric. A central question 
then concerns the role of the two orientations in the quantum theory. It 
has been claimed, for example, that the orientation of the triad reverses in a 
big-bang transition where the universe turns it inside out ([3], p. 67). This 
assumes that the universe is always in a state of definite orientation. But 
because of its linear structure, quantum cosmology allows the occurrence 
of arbitrary superpositions of the two orientations; one would then expect 
that the Universe is unlikely to be found in a definite orientation and that 
it is meaningless to talk about a change of orientation. How can this be 
understood? 

This situation reminds one of an old problem for chiral molecules posed by 
Hund in the 1920s [1]. Many molecules occur as objects with a definite shape, 
although the underlying Hamiltonian is parity-invariant. Energy eigenstates, 
which are superpositions of chiral states, are typically found only for small 
molecules, such as ammonia, but not for bigger molecules such as sugar. The 
reason for the occurrence of chiral states as robust quantities was understood 
in the last decades and can be traced back to the process of decoherence - 
the irreversible emergence of classical properties by the unavoidable interac- 
tion with the environment [SllSj- The spatial orientation of molecules is fixed 
by the scattering with photons and air molecules; the information about the 
superposition between different chiral states is transferred into the entangle- 
ment between the chiral molecules and the photons or air molecules and is 
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no longer accessible at the molecules themselves. 

In this paper, we shall examine whether decoherence is also responsible for 
the appearance of definite triad orientations in loop quantum cosmology. For 
this purpose, we need degrees of freedom that may serve as an appropriate 
'environment' for the triads, in analogy to the photons or air molecules for 
the chiral molecules. Such an environment must be of fermionic nature, 
because bosonic variables only interact with the metric and are unable to 
discriminate between the two different orientations of a triad. (We assume 
here that there is no parity violation in the matter coupling.) 

One possibility to implement fermions into quantum cosmology is to con- 
sider inhomogeneous fermionic fluctuations superimposed on a homogeneous 
background [7]. In fact, decoherence by fermions in quantum cosmology was 
discussed at length some time ago 0,(2]. It was found there that a fermionic 
environment leads to a suppression of interferences for the scale factor and a 
background scalar field and thus to their classical appearance, but in a way 
that is less efficient than for a bosonic environment; the reason for this smaller 
efficiency can be traced back to the Pauli principle. Inhomogeneous (bosonic) 
fluctuations can, of course, become themselves relevant as a quantum system 
subject to decoherence. This is the case for the primordial fluctuations in in- 
flationary cosmology, which can serve as the seeds for the origin of structure 
in the Universe; their decoherence is discussed in [10] and other papers. 

Here, we shall follow another route. Instead of taking inhomogeneous 
fermionic fluctuations, we retain homogeneity, but relax isotropy. Homo- 
geneity demands any flelds, including fermions, to be spatially constant. Im- 
posing, in addition, isotropy would have the consequence that for fermions 
this constant is zero, because otherwise one could construct from the spinor 
flelds a non-zero vector with a given direction, in contradiction to isotropy 
We therefore consider one of the simplest possible anisotropic models 
- the Bianchi I model (see, for example, [12] for a general introduction into 
classical and quantum aspects for such models). It is, in fact, the fermionic 
currents that in the considered model lead to the anisotropics. To simplify 
matters, we reduce this model further by imposing a rotational symmetry 
with respect to one axis, that is, by introducing the additional isotropy group 
U(l); one then arrives at a locally rotationally symmetric (LRS) model in 
which two of the diagonal components of the connection and the triad are 
equal. This was also done in [Tl]. The fermionic current aligned in the 
1-direction of internal space is then the sole reason for the anisotropy. 

The environment for the triads in the sense of decoherence is thus a 
homogeneous fermion fleld. It is known from quantum mechanics that one 
does not necessarily need many degrees of freedom for decoherence. For 
example, one can get decoherence from a long-range (e.g. gravitational) 
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interaction between two bodies, where one body decoheres the other (|13]. 
p. 228). 

In this paper, we shall basically use the model presented in [11] (fermions 
in a LRS Bianchi I model) and integrate over the fermions to study their 
decohering influence on the triads. But in order to investigate decoherence 
as a dynamical process, we have to introduce some notion of time. Quantum 
gravity is fundamentally timeless in the sense that there is no external time 
parameter present [1]. We thus introduce an inner variable and choose a 
massless scalar field v? for this purpose. Such a field was used as an inner 
time in loop quantum cosmology for the Friedmann case in [Hj and for the 
Bianchi I case (without fermions) in |15j . 

Section 2 contains a brief review of the formalism. Section 3 is the main 
part of our paper. We calculate the reduced density matrix and the linear 
entropy by tracing out the fermions and discuss the resulting decoherence. 
Section 4 contains our conclusions. 

2 The formal framework 

The formal framework of our discussion is based on the model of [11] sup- 
plemented by a massless scalar field. For the general formalism of fermions 
in loop quantum gravity, we refer to [TB] and the references therein. As for 
the units, we set c = 1 but keep G and h. 

In loop quantum gravity, the fundamental variables are constructed from 
the Ashtekar-Barbero connection A\{x) and the triad field Ef{x) [Il[3]. Here, 
z = 1,2,3 is the internal index of the 5u(2) algebra, and a = 1,2,3 is the 
usual space index. The fundamental variables are fluxes (integrals of the 
E°-{x) over two-dimensional surfaces) and holonomies (integrals of the Al^{x) 
along curves). For Friedmann-Lemaitre cosmology, the connection and the 
triad reduce to = c6l and Ef = pS"", respectively, where p = {a is the 
scale factor), and c is at the classical level given by c = 7a, where 7 denotes 
the Barbero-Immirzi parameterJl 

In the Bianchi I model, we have three scale factors instead of one. There- 
fore, instead of c we have three variables called Ci, 62, and C3 (we use here 
tildes to follow the convention of [H]); instead of p, we have p^, p"^, and p^. 
In the LRS Bianchi I model used here, we have 62 = C3 and p"^ = p^. 

The fermions are described by a densitized current J'a, a = 0, 1, 2, 3; we 
assume here, for simplicity, that = jTs = 0. The fermionic current is 

^The Barbero-Immirzi parameter is called /3 in [1], but we stick here to the convention 

of m]. 
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a source of torsion, which gives a new gravitational degree of freedom de- 
noted below by (j). The fermionic current depends on the four half-densitized 
Grassmann variables 0i, 02, 03, and 04 (summarized below under the label 
0). In the quantum theory, one has for the non- vanishing components of the 
fermionic current operator the expressions 

J^/h = 90,01 + 9e,02 - 90303 - 9e,04 (1) 

and 

Jr/h = 9e,0i + 9e,02 + 9e,03 + de,Q,. (2) 

In addition to the variables employed in [TT], we have here a massless 
scalar field, which serves as an inner time variable. The classical Hamiltonian 
for it reads 

H^ = MGpI{^\\p'\)-\ (3) 

(One recognizes that this Hamiltonian is insensitive to the triad orientation, 
which is why non-fermionic fields are not suitable to serve as a decohering 
environment.) For its quantization, we use the same technique as in 
The action of on a state is given by —ihd^. 

Since one of the fundamental variables in loop quantum gravity is the 
holonomy, the connection has to be exponentiated also in loop quantum 
cosmology in order to arrive at mathematically well defined expressions. (In 
homogeneous situations, the densitized triad components can be directly pro- 
moted to operators.) Functions on configuration space can then be expanded 
as follows: 

g{ci, C2, (p)= C'i.M2'*^ (^^^^1^1 + ^i/^2C2 + ik(pj , (4) 

where the sum is over finitely many /ii, /i2 G M, and G Z. It turns out that 
the Gauss constraint allows one to set /c = 0, so we shall omit k from now 
on. 

A general state \s) is defined in the full Hilbert space "H = "Hgrav ® 
"Hfermion ® "H^ and cau be expanded as 

1^) = 5Z ^'^I'/^^l©.'^) l/^l./^2) , (5) 

where /i2) denotes the common eigenstates of the triad operators and 
p^. Since fi2 ~^ ~/^2 corresponds to a triad rotation, we demand 

Smi,M2(0><^) = S;.l,-;.2(0,<^). (6) 
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In contrast to /X2, the sign of /xi determines the relative orientation of the 
triad. 

For the symmetric factor ordering, the difference equation for the total 
quantum state including the massless scalar field then read^ 



h{s^,,^,{e,ip)) :=2 [(|/i2 + 352| - |/i2 + 52|)Vl/ii + 25i| 

+ (1/^2 + ^2! - |/i2 - 52|)a/M S^i+25i,M2+252(07V5) 



+(1/^2 + ^2! - \fi2 - S2\)^A^J^1 

+2 



(|/i2 + 352I - |/i2 + ^2|)VlAti - 25i| 

S/ii-25i,M2+2<52(05 
(|/i2 - ^2! - |/i2 - 3(52|)V|Ail - 25i\ 

+ (1/^2 + ^2! - \fJ'2 - S2\)^A^^\ ■5^i-25i,M2-2<52(0,</5) 
-2 ((|/i2 - 52I - |/i2 - 352|)V|/il + 25i| 



+ (|/i2 + ^2! - |/i2 - 52|)Vl/il| ■S^,+25i,M2-252(07V5)) 



a/I/Wi + ^ll - a/Iah - ^1 

■ [(|/i2| + |/i2 + 452|)S/.i,/.2+452(0></') 
-4|/i2|S;,i,;,2(0><^) + (l/^2| + 1^2 -452|)s^i,/.2-4<52(0,<^)] 

(i/i2+52r/^-i/.2-52r/^)^ 



-^i^ir/^i^2r/^(i/ii+5ir/^ 



X 



4Y 



2Y 



1 + 7' 

1 + 72/12 



(3 + 272) 



7^ 



16 



;i2 



(7) 



The increments (5i and 62 appearing in this difference equation arise from 
the edge lengths of spin networks in the full theory; in minisuperspace, their 
exact form and meaning remains open. 



^For simplicity, we choose a — > cxd for the non-minimal coupling parameter a in 
which also leads to 6' = 1 and /? = 7 for the parameters 9 and /3 occurring in [TT]. This 
has no qualitative influence on our results. 
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We can write the difference equation ([7]) in the form 

27T ( J"^ J"^ f'^\ 
^(sm,M2(0,<^)) = ^ ( ci-^ + co-^ + c^-^J s^,,;,2(e,(^), (8) 

where we have introduced 

r := |/il|l/V2r/'(|/il + ^ll'/' - l/il - 5ir/')(|/i2 + 521'/' - |/i2 - 521'/')' (9) 



ci := --p-^ — ^, Co := ^ , c,^ := -16. (10) 



and 

1 + 72 ' ^° ■ 1 + 72 

We shall now discuss how superpositions of different triad orientations can 
be suppressed by the interaction with fermions. 



3 Decoherence of triad orientations 

In the following, we shall make the simplifications 62 = O3 = 64 = 0, 
01 = e and thus = and J^/tf = Jo/h = deQ = 1 - Qde- The full 
equation ([8]) can then be written in the form 

/(/^l, /i2)050s(/il, /i2, 6, if) = /i2, </:')s(/ii, /i2, 0, V^), (11) 

where the functions /(/ii,/i2) and (^(/ii, /i2, </?) encode the remaining depen- 
dences on /ii and /i2- We can then make the following ansatz for the solution: 

s(/ii, /i2, e, if) = So(Ail, /i2, ^) + 6si(/ii, /i2, if), (12) 

from which one gets the two equations: 

fsi = gsi, (13) 
gso = 0. (14) 

Since the total system is in a pure state, the total density matrix reads 

Ptot = s(/ii, /i2, 6', V5)s(/ii, /i2, e, V?). (15) 

Using the rules for Grassmann integration (see e.g. [IZ]), one can define the 
following reduced density matrix for the gravitational variables alone: 

Prcd(/ii, /i2; /"i, /U2; f)= dede e"®®s(/ii, /i2, e, f)s{fii, ^2, 0, </?). (16) 
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The integration yields 



Prcd(/Ul,/i2;/Ui,/U2;V5) = ■Sl(/ii,/i2,V5)Sl(/il,/i2,V5) + Sq (/i'l , /ia , Sq (/il , /i2, V^) • 



The total density matrix corresponds to an entangled state if both terms are 
present in the expression (JTTj) for the reduced density matrix, that is, if both 
So and Si are non-vanishing. Otherwise, the gravitational part is by itself in 
a pure state and there is no decoherence. 

A measure for the purity of the total state f|T5|) is the trace of p^g^j, which 
is equal to one for a pure state and smaller than one for a mixed state; it 
is directly related to the linear entropy Sn^ = 1 — trp^^^ [S] . One could also 
discuss the von Neumann entropy — A;Btr (prcd Inprcd)? but for the present 
purpose it is sufficient to restrict to 5*1111. 

To obtain concrete results, we calculate the evolution of the coefficients 
So and si in (|T^ . The spatial parameters pi and ^2 span a two-dimensional 
discrete lattice on which we have to choose appropriate initial data. Solving 
the equations in fll4p for the derivative in the continuous variable (f, we have 



Here, h{s) denotes the gravitational part of the difference equation, as given 
by the left-hand side of ([7]). 

To deal with (|T9|) . which is a combination of difference and differential 
equations, we use standard numerical techniques, implemented in a C source 
code written by ourselves. For this purpose, we have to discretize the con- 
tinuous variable ip. Since the equations in (IT^ are of second order in ip, 
we decompose them into a system of first-order equations. The resulting 
equations can be solved, for example, with the fourth order Runge-Kutta 
algorithm [18]. Now, we are in a position to calculate the coefficients sq and 
Si numerically, from which we then get immediately the functions s in the 
expansion of the state 

We now briefly describe the steps for the calculation of the reduced den- 
sity matrix. As initial data at <y9 = we have to specify sq, si, and their 
derivatives. We choose them to be Gaussians, normalized by the condition 
trpred = 1, except for Si{ip = 0), which is set equal to zero. This provides us 
with an unentangled state in the beginning of the evolution, see (fT7|l . Since 
Si vanishes initially, but not its derivative, the state will evolve into a mixed 
state. In the long-term evolution, it turns out that si will dominate over sq 
and the state will be pure again, as only the first term in f|T71) is left; this. 



(17) 





(19) 
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however, happens for times much longer than the times for which the model 
is applicable. 

The iteration at each time step starts with the calculation of sq and si, as 
described above. As a constraint on the numerical evolution, we normalize sq 
and si such that trpred = 1 is always preserved. Since the initial state is unen- 
tangled, trp^^^ is initially equal to one. As the inner time variable increases, 
the total state becomes entangled, and the purity factor decreases — the grav- 
itational variables are in a mixed state, and decoherence becomes more and 
more efficient. The result can be plotted as a function of the inner time 
variable ip, see Fig. [1] 
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Figure 1: The purity factor trp^gj plotted against the inner time variable ip. 

The interesting region is the right half of the diagram, where the purity 
factor decreases rapidly. In our simulation, the step size 6i resp. 62 is chosen 
in such a way that it is large in the beginning and smaller in the area of 
interest. As a test of the numerical robustness, we checked the algorithm 
with different numerical step sizes and different initial data. As long as the 
data are chosen as described above, the resulting curves look very similar. 
Another test is to choose the simpler Euler algorithm instead of the Runge- 
Kutta procedure. The qualitative behaviour is the same for both procedures. 

For the calculation, we have to fix the Barbero-Immirzi parameter 7 in 
([7]). In principle, 7 is a free parameter leading to a one-parameter ambiguity 
in the quantization. In most simulations, we have chosen 7 = 0.238, which 
is motivated by the desire to get the correct Bekenstein-Hawking entropy 
for black holes [19]. However, we have found that our results are largely 
independent of the specific choice for 7. 
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4 Conclusion 



In summary, one can say that fermions 'measure' the orientation of the triad 
in loop quantum cosmology. In this sense, they provide the universe with 
a definite orientation. The fermions act as an environment similarly to the 
photons or air molecules that can provide molecules with a definite shape. 
Fermions are thus an important ingredient in loop quantum cosmology. 

If the orientation is fixed by decoherence, it is hard to imagine a scenario 
in which the universe changes orientation by going from a pre- to a post- 
big bang scenario. One only has different branches of the total quantum 
state with different orientations. These branches are independent of each 
other, except possibly for small universes in which the decohering influence 
can be neglected, similarly to small molecules such as ammonia for which 
the influence of the environment is too weak to suppress the interferences 
between different chiral states. 

The situation may be compared to a scenario discussed in string quan- 
tum cosmology some time ago [2U]. In a semiclassical picture, the standard 
big bang is preceded in time (perhaps through a singularity) by a pre-big 
bang state. However, a consistent analysis in the quantum theory shows 
that wave packets cannot be continued through the singularity and that pre- 
and post- big bang just correspond to independent components of the total 
wave function that decohere from each other. Another example in quantum 
cosmology is the decoherence of the 'cosmological constant' (dark energy) by 
the interaction with metric perturbations [21] . 

The increase of entropy (decrease of the purity factor) arises because 
a special initial state is chosen. This initial state is characterized by the 
absence of entanglement between the gravitational and fermionic degrees of 
freedom. An initial unentangled state may be at the heart of the origin of 
irreversibility in our Universe [22l [23] . How this origin can be understood in 
full loop quantum cosmology is left for future investigations. 

Acknowledgement: We thank Martin Bojowald for helpful discussions. 
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